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Abstract 


A computer model has been developed which evaluates the performance of a heat 
exchanger. This model is general enough to be used to evaluate many heat exchanger 
geometry and a number of different operating conditions. The film approach is used to 
describe condensation in the presence of noncondensables. The model is also easily 
expanded to include other effects like fog formation or suction. 


Introduction 

It would be convenient in many situations to have a tool that could evaluate heat exchanger 
performance based on physical phenomena and general correlation only. This tool should 
stand alone, and should be able to predict heat exchanger performance in the absence of 
experimental data. This tool is needed as some heat exchangers (such as used in humidity 
control systems for space applications) might have little or no experimental data. Without 
this tool, prediction of performance for an exchanger becomes impossible. 

Existing commercial programs for calculating heat exchanger performance are usually too 
bulky. A source code of the program is not provided usually. Therefore, modifications of 
the code, which are often necessary to accommodate new effects, are next to impossible. 
This situation forced us to create a general program that can be customized easily for 
specific applications. In addition, the model can be included into simulation programs like 
Computer Aided System Engineering and Analysis (CASE/ A) or MATLAB widely used by 
NASA Ames scientific community. 

If there is no phase change taking place in a heat exchanger, the well-known Log-Mean 
Temperature Difference (LMTD) method is usually used to predict performance. For 
condensation of pure vapors, algebraic equations [1] are used to predict performance of the 
exchanger. For vapor condensation in the presence of noncondensable gases (air), a simple 
analytic solution is no longer possible. In order to describe energy transport in this 
situation, a reliable model must be developed. This model must be accurate and should be 
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easily expanded to include other effects such as fog formation or suction at the wall between 
the hot and cold stream to remove condensate. The model should also be easily customized 
to accommodate many types of geometry and flow conditions. 

Under some operational conditions condensation can start away from an inlet In this case a 
heat exchanger contains both dry and wet regions. The point between two regions is often 
described in the literature as a pinch point [2]. For dry region of the exchanger, the LMTD 
method can be used to predict performance. The LMTD method does not work for wet 
region in the case of condensation of vapor from vapor-air mixture. One of the most 
convenient way is to use a film model. 

The film approach provides a simple physical model for condensation of vapor-air mixtures 
which is widely used in engineering applications [3-5]. It was found that the film model was 
accurate in predicting heat exchanger performance in the condensation of steam-air mixtures 
[5]. In another study, the film model was used successfully to analyze condensation of 
binary mixtures [4]. The film approach has been also used in this model to determine 
performance characteristics for condensing vapor-air mixtures. 

During condensation of vapor from vapor-air mixtures, the noncondensable collects at the 
condensing surface as the vapor condenses [3]. This accumulation of noncondensable 
inhibits the flow of vapor to the interface and reduces the condensation rate. While the 
condensation rate is reduced, the heat transfer coefficient at the surface is actually increased. 
Diffusion induced velocity takes place as there is a higher concentration of vapor in the bulk 
than at the condensate film surface. Enhanced heat exchange was first predicted by 
Ackermann [4], An excellent review of past work in film model correction factors is 
presented in [4]. In [6], researchers observed an increase of heat transfer coefficient due to 
diffusion induced velocity experimentally. 
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A Description of the System Being Modeled 


The system being modeled is a counter-flow heat exchanger system, and a schematic can be 
seen in Figure 1. 

The hot stream consists of a vapor-air mixture while the cold stream consists of cooling 
water. The properties of the two inlet streams to the exchanger are known. For the air stream 
we know the inlet temperature, relative humidity, and flow rate. For the water side, we know 
the inlet temperature and flow rate. From this information, the performance of the exchanger 
and the properties of the outlet streams can be found. If there is no condensation in the 
exchanger the LMTD can be used to evaluate performance. If there is condensation in the 
presence of noncondensables, a more complex model must be used as LMTD theory is no 
longer accurate. 


Describing the System Past the Pinch Point 

In deriving the following equations, the common film approach assumptions have been 
made. The first assumption is that the hot stream has constant properties throughout its 
bulk. The second assumption is that all heat transfer takes place at the wall. In addition, the 
condensate film at the wall is assumed to be very thin. This is not a common assumption for 
the film model, but it allows the condensate and the wall to be at identical temperatures. The 
effects of vapor shear in the film layer are also disregarded. By applying the conservation 
laws of mass and energy to a differential area dA, the following equations are derived. 

AIR SIDE 

For temperature T h : 

C ph^^ = - h hVh-Tw) ( 1 ) 


3 



dfti dT 

Since — is many orders of magnitude smaller than — — , that term is ignored. Using 

dA dA 

this simplifying assumption, equation (1) is transformed to: 


Cph m h ~ hfi(Th T w ) 


( 2 ) 


Considering vapor concentration in air stream is negligible small, we assume that mass flow 
rate m/, is defined by dry air flow rate, which is constant We have then: 


For specific humidity H h : 



(3) 


W ATERSIDE 

For temperature T c : 

Cpctrtc ^- = -h c (T w -T c ) (4) 

WALL BOUNDARY CONDITION: 

h h (T h -T w )dA + H vap JdA = h c (T w - T c )dA (5) 


The boundary condition (5) is needed to balance the energy fluxes at the wall. In the 
boundary condition it is assumed that the wall provides no resistance to heat transfer. 

The air stream may reach a saturation at some point The bulk temperature T h can not go 

bellow a dew point temperature. Therefore, we assume that instead of (1) the bulk 
temperature T h is governed by an equation that followed from saturation curve. This 
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equation holds true downstream of the saturation point It corresponds to an assumption 
that the convective heat transfer removes from air an amount of heat that is required to keep 
the bulk temperature exactly equal to the dew point temperature. Thus, taking into account 
(3) we get equations for downstream of saturation point : 

For temperature T h : 

m h = —(Tsat)ff J (2’) 

dA 

where ( Tsat)n is the derivative of saturation temperature over specific humidity H. This 
function can be found elsewhere. 

WALL BOUNDARY CONDITION: 

C ph (Tsat) H JdA + H vap JdA = h c (T w - T c )dA (5’) 


Condensation Correlation 

To complete our system of equations, the rate of vapor condensation at the wall J must also 
be determined. Without the condensation term (J), the film model equations are reduced to 
an analytic solution from which LMTD theory is derived. According to the laws of 
diffusion, the condensation rate should be proportional to the vapor concentration difference 
between the bulk and the surface where condensation is taking place. In engineering 
practice, the Berman and Fuks correlation [7] has been widely used to predict condensation 
rates in the presence of noncondensables [8]: 

J = 0 . 52 -^ VRT(P WiA - P mp , w ) Vi P V \ 1 - P^apJ P)-°- 6 («) 
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where the coefficient of vapor diffusion D vap is related to the diffusion coefficient D by 
Dvap~ R S atTsat • using the Gibbs-Dalton ideal-gas mixture relations, and taking 
R sat fR h = 1.607 , equation (6) gives: 

J = ci n — VrT (SH h - SH w ) m {\ - \.miSH h y° 6 (7) 

UfyKsatl w 

where C\ =0.7 1 3 is the constant. However, the value C\ =2.495 is used for calculations, 
because this value provides the best fit of experimental data [17]. 

It was assumed in (7) that the vapor diffusion coefficient is calculated at a local saturation 
temperature, that is, at the wall temperature, because the vapor-air mixture at the wall is 
assumed to be saturated. Therefore the amount of vapor present in the mixture is then 
determined by the wall temperature [9]. Thus the rate of condensation at the wall is a 
function of the bulk specific humidity and the wall temperature. 


Heat Transfer Coefficient Correlations 


The program uses a Nusselt number correlation to estimate the sensible heat transfer 
coefficient The following correlations cover four different flow regimes on the hot and cold 
sides. 

For fully developed laminar flow, similar to [10]: 


,0677(Pr Re d hyi /L)' 

= Nun + — 

1 + 0.1 Pr (R ed hy( i/L) 


where Nuq is a theoretical constant which depends on a shape of tube cross section. For 
example, Nuq = 4.36 for a circular tube. For rectangular channels Nuq can be 
approximated by 

Nuq = 4.62(fl min /a max -l) 2 e -L63(a ™ i " / ^ ) + 3.61 
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where <2 m i n and <2 max are the minimum and maximum sides of the rectangular cross 
section respectively. 

For developing laminar flow [11]: 

Nu = 1.86(Pr Re 

For turbulent flow [10]: 

Nu = 0.0235(Re 0 ' 8 - 230)(1.8Pr 0 - 3 - 0.8)[l + (d hyd / L) 213 ] 

For entrance flow [15]: 

Nu = 0.33Re 0 ' 5 Pr 0J3 ( 

where the entrance length L ent is defined as 
^entj^hyd. ~ 0^5 ^ 


Miscellanies Effects 

Em efficiency: 

The program calculates also fin efficiency according to [12]: 
for radial fins with rectangular profile: 
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1 - 


2r re/ /i(arg)ATi(r rg ^ arg)- arg) Zx(r rg/ arg) 
(p(\ + r rel ) / 0 (r ret aig)K\ (arg) + tf 0 (r re/ ^g )A (arg) 


arg = 


0 

1 - r rel 


$ = (r e - r 0 ) 05 m 


m - 


2 ^ 
\ kS 0j 


1/2 


r re/ = 


_ r o 


where r Q and r e are radii of the base (outer tube diameter) and the tip of the fin respectively, 
k is the heat conductivity coefficient, 5 q is the fin thickness, and Ii( ), and Ki( ) (i=0,l) 
are i-Bessel functions. 


for longitudinal fin with rectangular profile: 
tanh mb 

v= ~^r 


where b is the fin height 
Pressure drop: 


Pressure drop is defined as: 


Ap=(f-L-+i.5)£L 


( 8 ) 


*hyd 


where / is the friction coefficient The value of 1.5 in (8) takes into account energy losses 
at an entrance and an exit of a tube. For turbulent flow (Re > 3000) the friction coefficient 
is [13] 
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0.221 


/ 


= 0.0032 + 7 rrz= 

Re 0 ' 237 


Capillary ri.K: 

As a result of a capillary force, condensed water may accumulate on the lower part of heat 
exchanger and block a portion of airpath. The blockage can be significant for heat 
exchangers with dense fin arrangement. The program calculates the height of capillary rise 
Heap • 



2(7 

Pgdfin 


where O' is the surface tension, djf n is the distance between fins, and g is the gravity. 

Note, that capillary rise formula is not applicable at microgravity. The capillary rise then is 
subtracted from fin height to determine actual vertical dimension of airpath. 


Cross-counter /low. 


Most heat exchangers employ cross-counter flow arrangement, which is less efficient than 
pure counter flow. The program determines the cross-counter flow factor T}^, which shows 

the reduction in heat transfer in a cross-counterflow heat exchanger vs. a pure counterflow 
heat exchanger [10] : 


n — 1 ^ 

% - 1 “ 7T 


r 

8 


with iff < 1 


_ Thin Jjwut 
N (AT m ) cont 


where N is the number of rows (in the air flow direction), (A T m 'f ont is the Mean Log 
Average Temperature Difference calculated for pure counterflow. Then, the heat transfer in 
the given cross-counterflow heat exchanger Q cross is; 
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qcwss = nkQ cont 


( 9 ) 


where Q cont is the calculated heat transfer for pure counterflow heat exchanger. 

7 *CTOSS r r , CVOSS 

toM and water T cout temperatures for cross-counterflow arft calculated 
then according to (9): 



( 10 ) 


where Tfoul and T^j are the outlet air and water temperatures respectively calculated 
for pure counterflow. For rare occasions when Iff > 1 , the program asks a user to use Fig. 
8-22 from [10] to determine the cross-counter flow factor T]^. and to correct outlet 
temperatures according to (1 1). 


Described above procedure is applied to a dry heat exchanger. In a case of condensation, an 
overall heat transfer coefficient multiplied by the cross-counter flow factor Tfy. is used in 

calculations. 


Program Method of Solution 

The flow chart given in Figure 2 describes the method of solution for the program. The 
region in the exchanger in which no condensation takes place at the interface is called the 
dry region of the heat exchanger. In the dry region of the exchanger, an analytic solution is 
used to determine the amount of heat transferred between the hot and cold streams [10]. In 
wet region of the exchanger (where condensation takes place), equations (2) - (5) are 
integrated to determine hot, cold side and wall temperatures and specific humidity profiles. 
Remind again, that the equations (2') and (5') are used instead of (2) and (5) for integration 
downstream of stream saturation point 
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The program integrates equations by using a Runge-Kutta integration scheme [14]. The 
Runge-Kutta method of integration solves an initial value problem, thus all conditions on 
one side of the exchanger must be known. As can be seen in Figure 1, on the one side of 
the exchanger the air inlet conditions are known and the water outlet conditions are 
unknown. In order to use the Runge-Kutta method, a water outlet temperature is estimated 
and the variables integrated down the length of the exchanger. The calculated inlet water 
temperature is then found and can be compared to the known inlet water temperature. It may 
take, more than one estimation to match the calculated cold stream inlet temperature to the 
actual temperature so the shooting method described in the flow chart is needed. 

Initially, the program calculates various constants that will be needed to find solutions in the 
dry portion of the heat exchanger (See Figure 1). These constants include the dew point of 
the hot stream at the inlet (Thdin), the heat transfer coefficients on the hot and cold sides, the 
overall heat transfer coefficient (U) for the dry portion of the heat exchanger, anticipated 
capillary rise, fin efficiency, and etc. Then the program determines if the entire heat 
exchanger is dry. This is accomplished by assuming that the wall temperature at the outlet 
of the exchanger (Twout) is at the dew point of the hot stream (Thdin). If the calculated cold 
stream inlet temperature (Tcinc) is less than the actual cold stream inlet temperature (Tcin), 
then the entire heat exchanger is dry and the analytical solution of [10] can be used to find 
outlet temperatures for the heat exchanger. 

If a wet section of the system is present, first the program checks to see if the entire 
exchanger is wet. The inlet wall temperature (Twin) is assumed to equal to the inlet dew 
point temperature (Thdin). Making these two assumptions, an inlet cold stream temperature 
is calculated (Tcinc) by integrating equations from the inlet to the outlet of the exchanger. If 
the calculated cold stream inlet temperature (Tcinc) is greater than the actual inlet 
temperature (Tcin), then the entire heat exchanger is wet The inlet wall temperature is then 
lowered until the calculated inlet temperature (Tcinc) is identical to the actual temperature 
(Tcin). 

If the exchanger is not entirely wet, then there is a pinch point somewhere in the exchanger. 
The program initially assigns the pinch point to the middle of the exchanger. The program 


11 



uses the analytic solution [10] to describe heat transfer in the dry section (air upstream of 
the pinch point). For the area in the exchanger where condensation takes place, equations (2 
or 2')-(5 or 5') are integrated to the end of the exchanger. Once the calculated cold stream 
inlet temperature (Tcinc) has been found, it is compared to the actual cold stream inlet 
temperature. If the calculated temperature is too high, the pinch point is moved towards the 
inlet If the calculated temperature is too low, the pinch point is moved towards the outlet 
This process is continued until the calculated inlet temperature (Tcinc) matches the actual 
inlet temperature (Tcin). 

The program then reports the outlet flow rate, relative humidity, and temperatures on both 
sides to the user. The program is designed to provide accurate results for heat exchangers in 
which the hot stream bulk temperature remains above or equal to its dew point. If the bulk 
temperature drops below its dew point, fog formation may start in the bulk [16]. This topic, 
and the adaptation of this program to include fog formation will be covered in a later paper. 


Comparison with Experimental Results 

The program has been used to calculate performances of a plate-fin heat exchanger made by 
Hamilton Standard (Windsor Locks, CN) [17] and a tube-fin heat exchanger made by 
Super Radiator Coils (Minneapolis, MN). The comparison of data is shown in Table 1. 
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Table 1 Experimental Data Compared to Theoretical Model 



"Hamilton 

plate-1 

Standard" 
fin HX 

"Super Radiator Coils" 
tube-fin HX 


Mrf. data 

Our 

results 


Our 

results 

Mrf. 
data ** 

Our 

results 

air flow (scfm) 

330 


224.0 


150 


air in temp. C° 

40.4 


19.0 


35.8 


air in r.h. 

% 

22.2 


67.0 




water flow 
Ipm 

7.62 


3.78 


1.89 


water in temp. 
CO 

6.4 


10.7 


26 


air out temp. 
CO 

11.3 

12.25 

11.8 

11.73 

26.8 

26.21 

air out r.h. 

% 


100 

100 

100 



water out temp. 
CO 

17.1 

16.82 

14.9 

14.83 

31.8 

31.83 

condensation 

rate/day 

23.8 

22.9 

6 

5.97 

0 

0 

sensible heat 
W 

5029 

4912 

931 

922 

766 


latent heat 
W 

683 

649 

180 

171 

0 


heat transfer 
W 

5723 

5561 

1111 

1093 

766 

769 


coil model 6 x 37.5 - 8R - 38/1 92 
coil model 4 x 50 - 7R - 38/1 44 


13 

































































Data in Table 1 shows that the developed program gives results very close to those given by 
manufacturers. It indicates that this general program can be used for different specific 
applications. 

The next improvement of the program is to be inclusion of fog formation. It allows us to 
determine more precise temperature distribution and condensation rate correspondingly. 
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Nomenclature 

A 


Appendix 


a 

b 

cl 

C f 

D 


dhyd 

D 0 

H 


h 


Hvap 

J 

k 

L 

m 

n 

MW 

P 

Pr 


Area [m^] 

Duct size 
Constant 
Constant 

Heat Capacity [fagfc] 
Diffusion Coefficient 
[ m /sec] 

Hydraulic Diameter [m] 

Pipe Diameter [m] 

Specific Humidity 

\kg(water)/ 1 

L /kg (dry air)J 

Heat transfer coefficient 

[%>*] 

Heat of vaporization [fc] 

Mass flux to wall \%/ 2 1 

L /m secj 

Heat conductivity coefficient 

[%d 

Exchange Length [m] 

Mass flow rate ^^sec] 

Constant 


Molecular weight 


yj 


Pressure [Pa] 
Prandtl Number 


R 

Gas Constant 






P 

Density 


Re 

Reynolds Number 

RH - 

Relative Humidity 

Sc 

Schmidt Number 

T 

Temperature [°C] 

U 

Overall Heat Transfer 


Coefficient 

\W/ . 1 
L Zm 2 K\ 

y 

Mole Fraction 



Subscripts used : 

c - Cold stream (water) 

d - Dew Point 

ent - Entrance 

g - Noncondensable gas 

h - Hot stream (gaseous) 

in - Inlet (actual) 

inc - Inlet (calculated) 

out - Outlet 

sat - Saturated 

vap - Vapor 

w - Wall 

Superscripts used : 

cont - Counterflow 

cross - Cross-Counterflow 
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Figure 1.- Heat Exchanger Schematic 
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Flgura 2.- Program Flow Chart 
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